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Abstract 



■ This paper concerns with the sensor management problem in collocated Multiple-Input Multiple-Output (MIMO) radars. After 
deriving the Cramer-Rao Lower Bound (CRLB) as a performance measure, the antenna allocation problem is formulated as a 

P-H ■ standard Semi-definite Programming (SDP) for the single-target case. In addition, for multiple unresolved target scenarios, a 

sampling-based algorithm is proposed to deal with the non-convexity of the cost function. Simulations confirm the superiority of 

C$ ' 

■ the localization results under the optimal structure. 

C/2 ■ 
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I. Introduction 
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o 

£ — ■ Multiple-Input Multiple-Output (MIMO) radars with collocated antennas have been introduced recently in the literature 

o; 

■ ifTSIITSl as an alternative to the traditional phased-array radar systems ll2D . Unlike the conventional phased-array radar systems 



in which the ttansmitted signals are fully-coherent, MIMO radars enjoy the orthogonality of transmitted waveforms. The 
orthogonality of the transmitted signals provides a number of benefits for the MIMO radars, such as the diversity in the paths 
0, virtual aperture extension [2 |, beam pattern improvement 0, and higher probability of detection [2|, over the phased-array 
systems. Consequently, this has generated much interest among researchers to analyze different aspects of collocated MIMO 
radars such as waveform selection [8||14|[17|[20|, range compression ifTTl and the applications of collocated MIMO radars in 
target detection, localization and tracking lfTTl lfT9ll25 1. 

The location CRLB was recently proposed in the literature as a tool to evaluate the localization performance of collocated 
MIMO radars [2||11||17|. The Direction-of- Arrival (DOA) of the target was defined as the parameter of the problem in (2). 
Then, the CRLB was derived according to the received complex signals. When multiple targets fall inside the same resolution 
cell of the MIMO radar, the CRLB might be also affected accordingly. An alternate form of the CRLB was then derived in 
IfTTl and the effect of the number of targets occupying the same cell might affect on the CRLB was analyzed. While previous 
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works only derived the CRLB for DOA estimation, it was shown in IfTTI that the range information of the target can be also 
included in the received measurements. Therefore, a novel measurement model was proposed in [11] and the CRLB was found 
for both range and DOA of the target. It was also shown that the CRLB is affected by the number and locations of targets 
falling inside the same resolution cell. 

Antenna allocation is a critical concern in MIMO array systems. An optimal antenna placement algorithm, where an array of 
closely-spaced antennas received the Time-of- Arrival (TOA) data, was proposed in Q. It is also shown in ||231 that the Posterior 
CRLB (PCRLB) E4l can be used to find the number and optimal locations of multiple sensors while there is no restriction 
on the closeness of inter-sensor distances. The CRLB was also employed in [9||13| for antenna placement in widely-separated 
MIMO radars. It was shown that the trace of the CRLB matrix can be written as a convex function of the location of antennas. 
Then, convex optimization techniques were applied in order to find the optimal placement of antennas. The CRLB was also 
used as a performance metric in iflOl for the antenna selection in widely-separated MIMO radars where a subset of antennas 
has to be chosen out of a large number of antennas that are widely-separated in the surveillance region. 

Recently, there has been interest in applying optimization techniques to different aspects of collocated MIMO radars as 
well. The CRLB was employed in ifTTll to find the optimal cross-correlation matrix of the transmitted signals where it was 
shown that the CRLB is a convex function of the cross-correlation matrix. A Gradiant-based approach was also formulated 
in (8) for beam-pattern synthesis by optimizing the transmitters' cross-correlation matrix. Although several other efforts have 
been made on the waveform optimization in collocated MIMO radars |5||6|, the problem of antenna allocation in collocated 
MIMO radars has not been addressed so far. In |[T3l , an algorithm for antenna selection in collocated MIMO radars was 
presented. Nevertheless, the proposed technique does not provide any systematic procedure for distributing the antennas in the 
surveillance region when the optimal set is chosen. It was demonstrated in fTTl that the CRLB of a collocated MIMO radar 
is a function of the location of antennas. The simulations in IfTTI also showed that the localization performance is affected by 
the distribution of antennas in the surveillance region. Consequently, it is of great interest to find an optimal distribution of 
antennas that provide the best localization performance. 

In this paper, the antenna allocation problem for collocated MIMO radar systems is addressed and a systematic approach is 
proposed based on the CRLB. To the best of our knowledge, there is no comprehensive work on the design and analysis of 
an optimal antenna placement framework for collocated MIMO radars. The main contributions of this paper are as follows: 

• A novel CRLB derivation for MIMO radars with collocated antennas: 

Although the CRLB was derived in the literature for the collocated MIMO radars [2 |[ 17], the effect of the range information 
was not considered in the CRLB derivation. In addition, there is no compact CRLB derivation in terms of the location of 



2 

antennas. In this paper, the CRLB is first derived for a collocated MIMO radar where both DOA and range of the target 
are both embedded in the signal model. Also, the impact of the situation in which multiple targets fall inside the same 
resolution cell is taken into consideration. 

• A convex optimization approach for the single-target case: 

It is shown that the antenna allocation problem can be dealt with by optimizing the location CRLB. To do this, the cost 
function is defined by applying suitable operators (e.g., determinant, trace, or maximum eigenvalue) to the CRLB. When a 
single target is located inside the resolution cell, the optimization algorithm is simplified to the well-known Semi-definite 
Programming (SDP) using the related convex relaxation techniques. 

• An optimization algorithm for the multiple unresolved target case: 

When multiple targets fall inside the same resolution cell, it is observed that the cost function is not convex anymore. In 
this case, due to the presence of sinusoidal terms in each entry of the Fisher-Information-Matrix (FIM), the cost function 
cannot be simplified into a convex form. Therefore, a sampling-based approach is proposed where initial conditions of 
the optimization algorithm are generated such that the algorithm moves towards the global minimum. Simulation results 
also confirm the efficacy of the proposed method in finding the optimum antenna allocation when multiple targets fall in 
the same or consecutive resolution cells. 

The rest of this paper is organized as follows. Section II presents a brief overview of MIMO radars with collocated antennas. 
The CRLB is derived for the MIMO system in Section III. Section IV deals with the antenna allocation problem where the 
convex optimization framework for the single-target case is described. Simulation results are given in Section V. The paper is 
concluded in Section V. 

A. Notations 

The notations used in this paper are as follows: 

• A = V(a): a diagonal matrix with An — and Aij = 0, i 7^ j 

• 3?(a): the real part of the complex variable a 

• 3(a): the imaginary part of the complex variable a 

• A/"(^i, £): a Gaussian function with mean \i and the covariance matrix E 

• T{A): the trace operator 

• A H : the Hermition transpose 

• A(:, i): the i-th column of matrix A 



II. MIMO Radars with collocated Antennas 
Consider an array of antennas with M transmitters and N receivers. 

Definition 1: Define Su = [xu yu]' and s r j = [x r j yrj] 1 as the location of the i-th transmitter and the j-th receiver in a 
2-dimensional surveillance region, respectively. 

Assumption 1: There are T targets in the region where x t = [x t yt]' denotes the location of the t-th target^j Also, the 
reflection of each target is modeled by a complex random variable a t = £ t + jCt with £ and ( being the real and imaginary 
parts of a, respectively. 

Assumption 2: It is assumed that the target's reflection follows a Swerling type I model [22 1 where {£ t ~ Af(£,t, a V)} an d 
{Ct ~ A/"(Ct, c^)}. Other models can be handled accordingly. 

Assumption 3: It is assumed that the distance between any two antennas is much smaller than the distance of the array to 
each target. It is also assumed that the arrays of transmitters and receivers are both collocated with the origin as the center of 
the mass of the array. 

Definition 2: Define h[fc] = [hi[k] ■ ■ ■ hM[k]] H as the transmitted waveform in the fc-th snapshot with K being the number 
of total snapshots. 

A. Signal Model 

Considering a collocated structure, resolution cells can be defined as a set of concentric circles where the radius of the c-th 
circle equals cr^n with run denoting the resolution width. Figure [T] shows a simple configuration of resolution cells as well 
as the antennas that are distributed uniformly. Note that the target is located inside the c-th cell. 

Assumption 4: It is assumed that all T targets are distributed in C consecutive cells (e.g., (c* + 1) to (c* + C)) where n c 
denotes the number of targets inside the c-th cell. Without loss of generality, it is assumed that c* = 0. 

Assumption 5: Transmitters send orthogonal signals with a diagonal cross-correlation matrix being defined as 

1 K 

R = -Y J mK H [k] = v{[p 1 ■■■ p M ]') (i) 

where P m denotes the total transmitted power by the m-th antenna. 

Definition 3: Defining rf = ||X(||2 as the Euclidean distance of the <-th target in the c-th cell to the origin, the ratio 
parameter /3 ( c is defined as follows: 

Pi = r * + (1 = C)nm (2) 

^ote that 3-D MIMO radars, although not very common in the literature, can be handled within our framework. 



4 




Figure 1. A collocated MIMO radar with three transmitters and three receivers. The target is located inside the c-th cell. The resolution cells are shown as 
concentric circles with crt,i n being the radius of the c-th cell. 

Now, given the above assumptions, the received output of the matched filter in the c-th resolution cell can be written as follows 
ED: 

E"=r(i-/3 t c+1 )# +1 c=o 

YZLxfttt c = C +™ (3) 

where w denotes a complex Gaussian noise with independent real and imaginary parts being distributed as {3t(iu), ~ 
A/"(0, o"^), and is the contribution of the t-th target in the signal received in the c-th cell, which is written as <\> c t = aftpf 
with the following form for the unknown term on the right-hand side of the equality ATI : 



Vc = < 



V>t = VKVBC(AzR?) (4) 

Here, WEC(A) stands for the matrix vectorization operator, and denotes the steering matrix of the t-th target defined as 
follows 07]: 

A\ = b c M) H (5) 
a t c = exp^-j^[sin(c9 t c ) cos(9 c t )}S t ) (6) 

/ 0_ s 



b t c = eX p(-j^[sm(0 t c ) cos(0 t c )]>S ) (7) 



where A is the wavelength, f?| denotes the DOA of the t-th target with respect to the origin, and the matrices St and S r are 
defined as 

S t = [sa ■ ■ ■ s tM ] (8) 
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Sri 



S r JV 



(9) 



Definition 4: Given the vector of the output of the matched-filter as 77 = [rj\ • • • 77^] H , define p = 
[Sfoi) 3(r?i) ••• Stfac) %(vc)}'. 

Now, the mean received output of the matched filter is defined as p = [8t(fji) Sfe) ' 1 ' 5ft(??c) ^(*7c)]'- The unknown terms 
5ft(?7 c ) and ^s(ij c ) can be found by calculating 5ft (0£ ) and 3(0^ ) as follows and then replacing in ([3J, respectively: 



(10) 



where the unknown terms on the right-hand side of the above equation can be written as follows: 



■ft(r;-| = \'Kr<,s ( y[sin(0 t c ) cos(6$)}Sl(S u S r ,R) 
9(^ c ) = VKsin ( — [sm(6 c t ) cos(8 c t )]n(S t , S r , R) 



(11) 



with O being defined as 



fl(5 t , S r , R) = (lixM 5 r -S t lixw) (iZ* 1i x jv) 



(12) 



where is the Kronecker product, and l a xb stands for a a x b matrix with all entries being equal to one. 

Given the signal model in (0 and the mean output of the matched-filter in ( fTOb , the following proposition provides the 
distribution of the output of the matched-filter ifTTl : 

Proposition 1: In a scenario with T targets located in C neighboring cells, the output of the matched-filter received by a 
collocated MIMO radar with M transmitters and N receivers (e.g., p) is Gaussian distributed with mean p and covariance £ 
defined as follows: 



/ 



S 



E11 


E12 





E21 


E22 






























\ 



J C(C-1) 



^cc 



(13) 



with the following definitions for E cc and E c ( c _ 1 ) terms: 

K<rl{Yr t Zx{P?? + °l)hMN c = C 
Kcrl{Y J n tUm i ?+ otherwise 

Yr t ;t\{i-P c t 2 ) 2 + °l)hMN 

Tl c 

^ c{c - 1) =KalY j {l-PtWthMN 



(14) 



(15) 



III. Cramer-Rao Lower Bound 

The CRLB provides the best Minimum Mean Squared Error (MMSE) bound for any unbiased estimator 1 1 1. In this section, 
the CRLB is derived for a collocated MIMO radar. It is shown that the CRLB is a function of the distances between any two 
antennas. Also, a scenario is considered with T targets distributed in C consecutive cells where different number of targets 
might be located inside each cell. 

Definition 5: For the £-th target located in the c-th resolution cell, define the state and parameter vector X t c and 9£, 
respectively, as follows: 



X 



(16) 
(17) 



The CRLB is the inverse of the Fisher-Information-Matrix (FIM) defined as follows (T): 

Definition 6: Assuming y as the received noisy measurements and 8 as the parameters of the measurement model, define 
the following matrix operator: 



Jee 



dlogp(y\8) f d\ogp(y\9)V 



06 



00 



(18) 



Refer to the definition of p and its distribution provided by Proposition 1 and define the stacked state and parameter vector of 
all targets as X = [(X{)' ■ ■ ■ (X*J (Xf)' ■ ■ ■ and 9 = [(9})' ■ ■ ■ (GjJ' (9?)' ■ ■ ■ respectively. In 

this case, the FIM is given by Jxx 1 - The defined FIM can be now written in the following form: 

d\ogp(p\X) ( d\ogp{p\X)\ 



Jxx' — E p 



dX 



dX 



(19) 



Using the chain-rule for partial derivatives, the above FIM can be simplified to the following form ifTTl : 



Jxx' = TJee'T' 



(20) 
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Here, T is called the system matrix and is written as 

7l 4 x4 
4 x4 72 

r = 



04x4 
04x4 



04x4 ' ' ' 4x 4 % 

with 04x4 as a 4 x 4 zero-matrix, and individual 7 t c terms being defined as 

|4 |4 

OX<7 OX? 



C 
nc 



(21) 



It = 



ae t 2K n n 

dyi dyt u u 



(22) 
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where the unknown partial derivatives can be derived using the definition of the ratio in (0 and the following equation for the 
target DOA: 

( vi 



— tan 



(23) 



Now, the FIM derivation becomes finding the unknown term Jee' m (f2Qb . The new Jee' can be broken into the following 
sub-matrices: 

•••0 



Jee' 



J(e 1 )(e 1 y ^(e 1 )^ 2 )' 



J, 



( 2)(Qiy J( 2)( 2y J(02)( e ay 



J, 







•^(e 3 )(e 2 )' >^(e 3 )(e 3 )' ^(e 3 )(e 4 )' 



(QC-i)(QC-2y J(0C-i)(0C-i)/ J(0C-i)(0C)/ 



- Jf, 







J 



(0 c )(0 c - l y 



J, 



(e°)(e°)' 



(24) 



Here, O c denotes a 4 x n c vector formed by stacking the parameters of those targets falling inside the c-th cell. The following 
equation can be written for 9 C : 

6 C = [(©D' ■•• (6')']' (25) 



Each individual entry in d24l can be also written as follows: 



^B'lVHfv ••• ^(e^xe^ 2 )' 



J (e? cl )(e^ 2 )' ••• J (e? ci )(e% 2 )' 

ci e {c 2 ,c 2 + 1} 



(26) 



x 



Finally, each entry of the FIM in (|26l ) can be simplified into the following form: 





fe^S 


"fesS 


*fec£? 






■fecS 




fe 1 C 2 


fe 1 0% 


■fees 


fees 


"fe^S 


"fe/3^ 2 


fees 


■feeS 



(27) 



Note that the matrix given by ((27) is a 4 x 4 FIM sub-matrix that includes the information correlation between the parameters 
of the n-th target in the ci-th cell and the m-th target in the C2-th cell. Before presenting the algebraic expressions for each 
entry of the FIM given in ( f27l i. the following new notations are defined: 

Definition 7: Assuming that the n-th target is located in the ci-th cell, define the following new notations: 



<(0 = yNK 1 ) cos(0£)]fiO,O 
= [cos(^) -sin(^)]' 

where fi(:, Z) denotes the Z-th column of matrix ft with Sl(St, S r , R) being written as ft for brevity. 

Definition 8: For any two targets falling inside cells c\ and c%, respectively, the following notations are defined: 



(28) 
(29) 
(30) 



(31) 
(32) 



nm _ rtiAC 2 _/-ci£C 2 
c±C2 >n Sm Sro 

The covariance matrix £ found in ( fT3l can be now rewritten for cells {ci — 1, c\, c 2 } with cx G {c2, C2 — 1}. Using the general 
form given by ( fT3l and expressions provided by ( fT4b . the new covariance matrix can be written as 

ci C4 



c 4 c 2 c 5 © h\iN 
c 5 c 3 

where c% terms are found using (TT4l . Similarly, the new notation is defined as 

P* - ^fe-i) RfoJ 3(r? Cl ) Kfe 2 ) 3(t7c 2 )]' 

Now, it can be shown that the inverse of S* can be written in the following form: 

k\ fc 4 k$ 
fc 4 fc 2 fc 6 
fc 5 fc 6 fc 3 



(33) 



(34) 



£7 X = 



I2MN 



(35) 
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The following proposition provides algebraic expressions for each entry of the FIM in d27| i: 

Proposition 2: Assume a scenario with T targets falling inside C consecutive resolution cells. Each entry of the FIM defined 
by (|27T i can be calculated as follows: 



2n 



MN 



53(jff)'n(:,i)fi'MK 



i=i 



x x co« (I) - < (/)) + t ™ sh« (I) - u# (/))) 



C fl <=i, 



bin 



11=1 



Jgiffl = K 



Cn 



K 2vr 

r bin A 



MN 



K 



Jfl c l? c 2 = it" 



2-k 

\ 

2tt 



K 
K 

i bin 

K 



53 <™ cos« (0 - u$ (/)) + c 2 sin « (0 - < (0) 

UN 

^cos«(0-<(0) 
(=1 

53 (P^)'n(:, { <? 2 cos« 2 (0 - < (0) + 6- sin« (J) - w « (/))} 
i=i 

MN 

5>?)'«0> {-C? co«(z) - w * (/)) + e sm(^ (o - w « (0)} 

(=1 

MTV 

5>°W,0 {ecos«(0 -w«(0) + C sinK?(0 -<(0)} 
{=1 

fAT 

53 {C 1 cos«=(/) - + C 1 sin«(/) - <(!))} 

^ {C 1 cos«(0 - <(0) - e sinK?(0 - <(!))} 

IV 

53 sin«i (l)-<(0) 



A 



Ll=l 

MN 



Cpl 1 el 



Cpl 1 el 



i=i 



(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 



with F(/3!£, (3%l) being a known function of ratios, and the following expression being given for unknown coefficients in the 
right-hand side of the above equations: 



C p% Pi 



'(My o" 


z- 1 


fe'" 




K 1 


o fey' 



Cl = C2 - 1 

otherwise 

[-1 1 0] S- 1 [0 -11]' ci = pa - 1 
[0 -11] Ej 1 [0 - 1 1]' otherwise 



3^ ' °/3^6l 



(45) 



(46) 



(47) 



0#2 ' "» 

Proof: See Appendix lAl ■ 
Note that the above proposition can be used to find the FIM for every pair {O^ 1 , ©„} where 1 < {ci , C2} < C and 1 < n < n Cl 
and 1 < m < n C2 . The FIMs calculated in ( T2"7l i are then inserted in ( T26b and d24"l i. respectively, to obtain Jee'- The CRLB is 
finally found by inverting the FIM as 

Cxx> = (T-^'Cee'T- 1 (48) 
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IV. Optimal Antenna Allocation 

A. Motivation 

It can be shown that the localization performance of the collocated MIMO radar is affected by the distribution of antennas 
in the surveillance region. To illustrate this, consider a representative scenario with two antennas (N = 2, M = 2), where each 
antenna can both transmit and receive signals. We take a single target scenario into consideration with parameters [30° .33 1 1]', 
which is located in {r, 6} = [825m 30°]'. The variance of DOA estimates {Cg2) is now shown in Figure |2]in terms of different 
inter-antenna distances for the designed scenario. It can be observed that the geometry of sensors (inter-sensor distances) 
affects the performance bound of DOA estimation, where the estimation variance at the minimum point is 33% lower than the 
maximum variance. 

Unfortunately, the graphical tool cannot be developed for cases with more antennas. Therefore, this section concerns with 
designing a systematic algorithm for the antenna allocation problem in collocated MIMO radars. First, the case with a single 
target in the surveillance region is considered. It is shown that by considering suitable geometric constraints, the antenna 
allocation problem can be formulated as an SDP procedure [4 |. Then, the problem is extended to the case with multiple targets 
in the same or consecutive resolution cells. It is shown that the derived cost function is non-convex and a sampling-based 
approach is proposed to capture the global minimum of the cost function. 
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B. Single Target Case 

When a single target is placed in an arbitrary resolution cell, all cos(.) and sin(.) terms in the individual entries of the FIM 
defined by d27b vanish. Let us assume 8 = [9 C j3 c £ c £ c ] ' as the parameter vector of the single target fallen in the c-th cell. 
Using the results given in Proposition [2] it can be observed that only terms Jqc^c are a function of the antenna locations where 
v c € {9 C , (3 C , £ c , £ c }. On the other hand, according to the definition of the matrix SI in 0121 . one can show that: 

MN 



5^O(:,0 = (49) 



It can be observed that only the term J^c-p can be considered as a function of the antenna locations. 

Definition 9: Define the difference between the m-th transmitter and the n-th receiver as As nm = St m — s rn . 

Corollary 1: In a collocated MIMO radar with M transmitters and N receivers, where a single target is located in the c-th 
resolution cell, the FIM is a function of inter-antenna differences. In addition, all entries of the FIM are independent of the 
inter-sensor differences except Jrgc\2, which is also a convex function of the unknown differences. 

Proof: It was shown that only J{e c ) 2 is a function of the sensor locations. Now, it is demonstrated that it is a convex 
function of the parameters (difference vectors). Using the algebraic terms given by Proposition 2, the entry J(ec)2 can be 
simplified into the following form: 



■In = l<\ ^) \a c \ 2 



A 



5> C )'0(:,Z)0'(:,OP C 



C ( 0cy (50) 



.i=i 

Consider the definition of VL in (fT2l . It can be then observed that f2(:,Z) is a linear function of the corresponding difference 
vector As„ m . It is also known that J^g^2 is a convex function of 0(: , Z) terms due to the appearance of quadratic terms in 
(l50l [4|. Therefore, J(gc)2 is also a convex function of the difference vectors. ■ 
The antenna allocation problem can be now dealt with by minimizing the trace of CRLB, maximizing the determinant of FIM, 
or minimizing the maximum eigenvalue of CRLB ifTTl . The following lemma proposes the convex optimization formulation 
for the antenna allocation problem in a collocated MIMO radar system where a single-target scenario is considered: 

Lemma 1: Consider a collocated MIMO radar with M transmitters and N receivers. In addition, assume that there is a 
single target located in the c-th resolution cell. Then, a convex optimization algorithm that finds an optimal placement of 
antennas is given as follows: 

max(A Sll ,.,AS„ m ) J{0"Y (51) 

Proof: The optimization problem can be formulated as minimizing the determinant of the CRLB, which is equivalent to 
maximizing | Jxx'\- In addition, the system matrix T defined in (|2TT > is independent of the location of the antennas. Therefore, 
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the final goal is to maximize |J©o/qo)/|. Now, the FIM in d27l ) can be written in the following new form: 



J(gc)2 b' 



B 



(52) 



where b and B are block vector and matrix, respectively, formed by remaining entries of Je c e c in ( |27| >, respectively. The 
determinant term can be written as 

|J e c 0c | = |fl||J (fl . )a -b'flb| (53) 

It is known that both B and b are independent of the antenna placement. Therefore, the determinant maximization can be 
achieved by maximizing J^y with respect to ft. However, it is also known that O is a linear function of As nm terms. The 
optimization problem can be finally simplified to maximizing J(gc)2 with respect to the As„ m terms, which is the final form 
given in (IBTl l. ■ 
The final optimization problem can be now constructed by imposing the following constraints on the inter-antenna distances: 
The inter-antenna distance: 

In practice, antennas need to be well-separated to ensure maintenance and safety considerations. In addition, the inter-antenna 
distance should be small enough to have the far-field assumption still valid. Based on the given targets, the following constraints 
can be considered: 

||As nm || 2 > d nm (54) 
||As„ m || 2 < e nm ,V m = 1,...,M, n = 1,...,N (55) 

where e nm and d nm are design parameters. 
The center of the mass constraint: 

It was mentioned in Assumption [3] that the center of the mass of the array is located in the origin. Therefore, the following 
new constraints are formed on the location of antennas: 

M N 
m— 1 n— 1 

Note that the FIM is a function of inter-antenna distances and therefore, a set of optimal difference vectors might correspond 
to an infinite number of sensor locations. The constraint given by d56i l ensures that the mass center of the obtained geometry 
is in the origin. The uniqueness of optimal solution is further discussed in this section. 
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(57) 



Considering the above defined constraints, the new optimization problem can be written as follows: 

maX {ASlli ... )ASnm} £m=l E^l(P C )' As «™ As nmP C 

S.T ||As nm ||2 > dnm 

||^S nm ||2 ^ &nm 

Em=i s *™ + En=i s ™ = 0, V m = {1, • • • , M}, n = {1, • • • , iV} 
In writing the above equation, it is assumed that the transmitted powers are all the same and unitary (Pi = P 2 = • • ■ = Pm = 1). 

The optimization problem given by d57l i is not convex and therefore cannot be solved using the standard approaches. The 

following theorem reformulates the optimization problem in ( l57l i as an SDP: 

Theorem 1: Consider a single-target scenario with a collocated MIMO radar being used as the measurement tool. Defining 

T* = {Tn, • • • , T nm }, S* — {sti, • ■ ■ , s t n}, and t — [tn ■ ■ ■ t nm ]', the optimal placement of transmitters and receivers that 

maximizes the determinant of FIM is found by solving the following SDP optimization problem: 



maxj,, s , j- 
S.T. 



EM v^-^ + 
m=l L^n=\ Tnm 

Em • i -^ N _ n 

m=l S *m + L~in=\ S "i — U 



P (PnmP) ^ t Tl 



I2/.2 Strn S rn 
s' - s' -P 2 



^0, 



Si 



3 tm a rn 



< 



(58) 



c' _ c' 
"tin a rn 



^0, Vro = {l,---,M},n = {l,---,tf} 



Sim S rn Pnm 

with P = p c (p c )', and ^ as the generalized inequality operator. 

Proof: See Appendix IE1 ■ 
The above optimization problem can be now efficiently solved using standard packages for solving SDP problems fl2l . 

Remark 1: The optimization problem in ( fSTt proves the dependency on the parameters of the target through the matrix P. 
The following proposition shows how the optimal structure is affected by changing the DOA of the target: 

Proposition 3: Consider a single-target scenario with a collocated MIMO radar being used as the measurement tool. Defining 
0i and 62 as two different DOAs and {S° 1 ,S° 1 },{S° 2 ,S° 2 } as the assigned optimal antenna allocations, respectively, the 
following equations are valid: 



„o2 
'tm 



„<>2 



GA0S° m 



G Ae Kn V m = {1, • • • , M}, n = {!,■■■ , N} 



(59) 
(60) 
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with A6 = 02 — 61 and Gab as the rotation matrix defined as follows: 

cos(A(9) -sin(A0) 

Gab 



(61) 



sin(A0) cos(A6>) 

Proof: See Appendix ICl ■ 
Remark 2: The SDP formulation given by Theorem 1 does not provide any information regarding the uniqueness of the 

optimal solutions for the location of antennas. The uniqueness of solutions is now discussed in the following proposition. 
Proposition 4: Consider a single-target scenario with a collocated MIMO radar being used as the measurement tool. Then, 

there are at least two solutions for the optimization problem in ( TSTT i as {S° , S® 1 } and {S° 2 , S° 2 }. In addition, the first optimal 

configuration can be obtained from the second one by a simple rotation as follows: 

s t °? = G^°l i = {l,---,M} (62) 
s°| = G^ 2 , j = {l,---,N} (63) 

where G^ is a rotation matrix with ir as the angle of rotation. 

Proof: See Appendix ID1 ■ 

C. Multiple-target Case 

When multiple targets fall inside the same resolution cell (or consecutive cells), the individual entries of the FIM in d27b 
are no longer convex. 

Proposition 5: Consider a collocated MIMO radar system with M transmitters and TV receivers with d nm < ||As„ m || < 
e nm , V m = {1, • • • , M}, n = {1, • • • , N}. Also, assume a scenario with two targets in the c-th resolution cell with parameters 
Of and G>2, respectively. Then, the term (ujf(:,l) — £a>2(:,Z)) falls in the following interval: 



d nm ^2{\ - cos(^ - 91)) < I) - uj c 2 (:, I) < -^e nm ^2 (1 - cos(0 2 c - 0f)) (64) 



with 1 = {1,---,MN}. 

The above proposition states that the more separated the DOA of targets, the wider I) — uj^i', I))- For example, defining 

Aui as the difference between the upper and lower bounds of I) — uJ^i'-i 0) i n < l64l >. Figure [3] shows how Aui changes by 
varying the difference between the DOA of targets. It is observed that when the targets are well-separated in the DOA space, 
the difference between the maximum and minimum bound is significant. This also highlights the contribution of the sinusoid 
terms in each entry of the FIM, which might result in several local optimum points. On the other hand, the convex relaxation 
approach used for the single target case cannot be applied to the cost function derived for the case with multiple targets in the 
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Figure 3. Variation of the sinusoid argument by changing the difference between the DOA of targets. 

same cell. The above problems make the optimization problem non-convex when there are more than one target inside each 
resolution cell (consecutive cells). 

To handle the above problem, the optimization algorithm is solved for different initial locations of the antennas. However, 
a large number of initial points are required to capture the non-convexity of the cost function. The sampling approach is now 
proposed in Algorithm Q] In the proposed algorithm, Q denotes the covariance of the normal density function that is used to 
generate new initial points. While the covariance matrix is chosen experimentally, a small variance might make the algorithm 
be trapped in the local optimum point. Therefore, an intelligent choice of the covariance matrix can enhance the efficiency of 
the algorithm. The main idea behind the proposed approach is to, first, find an estimate of the optimal antenna location, which 
might be a local solution. Then, initial points are generated based on the obtained optimal location. 

Remark 3: Note that the proposed algorithm does not always guarantee that the optimization algorithm captures the global 
solution. Due to the non-convexity of the cost function, there is also no analytical way to capture the global solution. 
Nevertheless, the proposed sampling approach initializes new points around the initial local solutions and pushes the overall 
algorithm into the global solution. As shown in the simulations, it can be observed that regardless of the initial selection of 
the antenna location, the algorithm always converges to a unique solution. 

Remark 4: The procedure given in Algorithm [T] is terminated when the cost function is not reduced in more than /i iterations 
where the parameter [i is empirically chosen. If the algorithm finds the global solution of the cost function, randomly-generated 
initial conditions around the optimal point does not give a lower cost and, therefore, the algorithm does not advance in subsequent 
iterations. In this case, the procedure is stopped after \i unsuccessful trials. 
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Algorithm 1 The optimization algorithm for the case with multiple targets in the same resolution cell 

Initialization: Generate an initial location of antennas as s° m , S° n with m = {1, • • • , M} and n = {1, • ■ ■ , N}. 

Optimization: Find an optimal distribution of antennas by minimizing the following cost function: 

min Stl ,... : s rJ v T(C X x>) 
S.T. ||As nm || > d nm 

||As nm || ^ &nm 

Eti s te + E!=iSrn = 0, Vm = {l,--.,M}, » = {1,...,JV} 

Optimal Cost: Initialize s° m and s°„ and calculate the assigned cost as C° = T(Cxx , )s° 1 ,—,s° N - 
Sampling: While u < U or NA < /i: 

. Sample s° m - Af(s% m ,Q) and s°„ - A/"(s°„, Q) with m = {1, • • • , M} and n = {1, • • • ,N}. 

• Run the optimization algorithm and find the new distribution of antennas S( m ,s*„ and associated cost C*. 

• if §tt < 1 then 

s°tm = Sf TO , s° rn = s*„ with m = {1, • • • ,M} and n = {1, • • • , N}. 
C° = C*. 
NA = 0. 
. else 

NA = NA + 1. 

• end if 

Report s° m and s° n as the optimal distribution of antennas. 



(65) 



V. Simulation Results 

In this section, we analyze how the optimal allocation of antennas in the surveillance region affects the localization 
performance of the MIMO radar system. To do this, a collocated MIMO radar is first designed with the parameters in Table U 
In the following subsections, the performance of the optimization algorithm is studied first for a single target scenario. Then, 
the simulations results will be provided for a scenario with multiple targets occupying the same resolution cell. 

A. A Single-Target Scenario 

Initially, consider a single target located at [410 — 710]' (m). The parameters of the target are also chosen as follows: 

6 = [~ .33 3 3] (66) 
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Table I 

Simulation Parameters 



Parameter 


Description 


Value 


Tmax 


Maximum coverage range of transmitters 


5 (km) 




Range width 


30 (m) 


A 


Wave-length 


30 (cm) 


K 


Number of snapshots 


128 




Variance of the scatterers 


io- 4 


°i 


Variance of the additive noise 


1 


Pm 


Transmitted power 


1 (W) 



In the first experiment, assume that there are M antennas where each antenna can both transmit and receive signals. Two 
antenna configurations, a Uniform-Linear-Array (ULA) with half wave-length spacing and the optimal geometry proposed 
in this paper, are considered in this part. For simulations, it is assumed that d mn = A, e mn = 2A V{m, n}. The optimal 
configuration of antennas is shown in Figure |4] for different number of antennas. In addition, Figure [5] presents the location 
CRLB for both optimal and ULA structures separately. It can be observed that the CRLB of the optimal configuration is much 
lower than that of the ULA structure. The improvement becomes more significant when the number of antennas is smaller. For 
example, for the case with M — 2 antennas, the CRLB of the optimal structure is around 6 times lower than that of the ULA 
configuration while the improvement decays to 2 times lower at M — 5 antennas. When the number of antennas increases, 
the gap between the optimal and ULA CRLB becomes smaller because the Signal-to-Noise Ratio (SNR) is large enough to 
make up the poor geometry of antennas. 

1) The effect of the target DOA on the optimal structure: Consider the above scenario with M = 4 antennas. While the 
target is still assumed to fall in the same cell defined in the above experiment, its DOA varies in the interval [-^-, §]• The 
optimization algorithm is now implemented to find the optimal configuration of antennas. Figure [6] shows the results for four 
different target DOAs. The results shown in Figure [6] imply that the optimal configuration with 6\ as the DOA can be obtained 
from the optimal structure with 62 by rotating the geometry (62 — 6\) (rad) around the mass center, which confirms Proposition 

13 

2) The localization performance of the optimal structure: Assume M = 3 for the number of antennas. Besides the optimal 
and ULA configurations, a random antenna allocation is also used for the test where the antennas are randomly distributed in 
the underlying surveillance region. The localization Root Mean Squared Error (RMSE) is now calculated at different target 
SNRs where all results are obtained after 100 Monte Carlo simulations. Figure [7] presents the resulting RMSE for each of 
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Figure 4. The optimal configuration of antennas for a single-target case. The optimal configuration is found for different number of antennas where each 
antenna can both transmit and receive signals. 




Figure 5. Localization CRLB for the ULA configuration and the optimal structure. The CRLB is found for a single-target scenario and different number of 
antennas. 

the above configurations. It is observed that the optimal configuration achieves the lowest RMSE while the ULA provides the 
worst results. The random allocation also gives an RMSE between the optimal and ULA configurations although other random 
distributions of antennas may provide higher RMSE results. 

3) The optimal design for separate transmitter and receiver arrays: Simulation results are now provided for a scenario 
in which each antenna can either transmit or receive signals. Consider a single-target scenario with 6 C = — j (rad) as the 
DOA. The optimal structure is now found for two cases with M = N = 2 and M = N = 6 antennas. Figure [8] presents the 
obtained optimal structures where, for each case, the results are given for scenarios with the same and separate transmitters 
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Figure 6. The optimal configuration of antennas for a single-target case, and for four different target DOAs. 

and receivers, respectively. It can be observed that the optimal structure obtained for each case (e.g., the same and separate 
transmitters and receivers) is the same with transmitters and receivers being clustered such that the mutual distances between 
the same-type antennas (e.g., transmitter or receiver) is minimized. To test this hypothesis, assume that 6 antennas are available 
and there are two scenarios with M = 4 and M — 3 as the number of transmitters at each scenario. The optimal structure is 
now found for each scenario and the final results are shown in Figure [9] It can be observed that the same optimal structure is 
obtained for both cases with antennas being clustered based on the mutual distances between the antennas with the same type. 

Note that although the obtained optimal structures in Figure [9] are similar, the optimum cost function might be different 
based on the number of signal paths (M x N). For example, for the configurations given in Figure |9l the optimal cost is 
calculated to be 0.7545 and 0.6393 for M = 4 (M X N = 8) and M = 3 (M X N = 9) antennas, respectively. The obtained 
optimum cost values also confirm the fact that the more the diversity gain, the lower the achieved optimum cost. 



B. Multiple-Target Case 

In this subsection, the optimization algorithm is applied to a scenario with more than one target being located in the same 
resolution cell. Let us assume there are two targets falling in the same cell with the following parameters: 



&l = 



e c 2 = 



-- .33 3 3 



.66 3 3 



(67) 



Based on the results in Figure [3] it is now evident that the effect of sinusoidal terms on the cost function cannot be ignored 
due to the large value for A6. First, the optimization framework given by Algorithm Q] is applied to the two-target scenario with 
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Figure 7. Location RMSE for different target SNRs and for the MIMO radar with M = 3 antennas. The RMSE results are obtained for three different 
structures (ULA, optimal, and randomly-distributed configurations). 
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Figure 8. The optimal configuration of antennas for a single-target case. The optimal configuration is found for different number of antennas where each 
antenna can either transmit or receive signals (Q- transmitters, □- receivers, and A- transceivers). 



different initial conditions. Figure \TU\ shows the cost values obtained at different iterations of the algorithm and for different 
initial conditions. It is observed that the algorithm captures the global minimum after a number of iterations. While each initial 
condition leads to a different cost value, the sampling approach finally finds the structure corresponding to the global minimum. 
Note that without the sampling procedure, each initial condition leads to a different optimal cost as shown in Figure QT| 
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Figure 9. The optimal antenna configuration for the single-target scenario with M + TV = 6 antennas. The optimal structure is found for two cases with 
M = 4 and M = 3 antennas as the number of transmitters (Q- transmitters and □- receivers). 



Cost curves for different initial conditions 




3 4 5 6 

Optimization iteration 



Figure 10. The cost for 10 different initial antenna locations (different symbols correspond to initial conditions). The simulations are done for a two-target 
scenario the same resolution cell. 

1 ) The effect of the angular separation on the optimal configuration: The optimal structure is found for different values of 
AO = #2 — 0\. The optimal configurations are now depicted in Figure Q~2] for four different values of A8. It is observed that 
when AO — > 0, the obtained structure resembles the one given in Figure [4] for the scenario with M = 4 antennas. Nevertheless, 
for other values of AO, a new structure whose geometry depends on the distribution of targets in the resolution cell is obtained. 
Figure Q~3] also presents the cost function (i.e., the trace of the location CRLB) for different values of AO where the results are 
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Figure 11. The optimal cost for 10 different initial antenna locations without using the sampling-based approach. 

obtained for both optimal and ULA configurations. The graph indicates that the closer the targets the poorer the performance. 
For example, the cost at A9 — is 100% higher than the one at A9 = when the optimal structure is taken. In addition, 
Figure Q~3] confirms the superiority of the optimal structure to the ULA. The rate of the improvement also increases when 
targets become more closer, with 10 times lower cost at A9 = ^ compared to 5 times lower cost achieved at A9 = 2^. 

The obtained results in Figures [4] and Q~3] imply that although the optimization algorithm can be implemented more efficiently 
when the angular separation between two targets becomes smaller, the performance of the localization is degraded. In other 
words, there is a trade-off between the quality of the localization and the efficiency of the optimization algorithm. Smaller 
values of A9 makes the FIM entries in d27i > less dependent on the sinusoid terms. 

2) The optimal design for multiple unresolved targets: The performance of the optimization algorithm can be also evaluated 
for a scenario with more than two targets inside the same resolution cell. It is known that there is a bound on the maximum 
number of targets that can be uniquely detected in the same resolution cell fl6l . Assume different number of targets are 
placed in the c-th resolution cell with the same SNR being assigned to each target. Also, consider the MIMO structure with 
M = N = 4 antennas where each antenna can both transmit and receive signals. We find the optimal structure for each case 
with a different number of targets inside the same resolution cell. 

For comparison, the localization algorithm is also applied to the obtained structures and the location RMSE is calculated 
by taking an average of individual estimates in 100 Monte Carlo runs. The RMSE results as well as the location CRLB 
are now depicted in Figure [14] where the graphs for the case with the ULA MIMO structure are also included. While the 
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Figure 12. The optimal antenna configuration for the two-target scenario. The optimal structure is found for different values of AO. 

localization performance degrades, when the number of targets increases, the optimal structure always shows the lower RMSE 
compared to the ULA configuration. In addition, when more targets fall inside the same resolution cell, the difference between 
the obtained RMSE of the ULA structure and that of the optimal configuration becomes higher. For example, for the scenario 
with T = 2 unresolved targets, the optimal RMSE is 53% lower than the RMSE obtained by the ULA structure. Nevertheless, 
the gap widens to 123% when 5 targets occupy the same resolution cell. Although the distribution of targets in the cell also 
affects the localization performance ifTTI . this experiment shows the superiority of the optimal structure compared to the ULA 
configuration, specially, when more targets are placed in the same resolution cell. 



VI. Conclusions 



This paper considered the antenna allocation problem in a collocated MIMO radar system. A novel derivation of the CRLB 
was presented where both range and DOA information were included in the CRLB. An SDP problem was then formulated for 
antenna allocation when a single target is located inside the resolution cell. Then, the antenna allocation was extended to the 
multiple unresolved target scenarios, and it was shown that the final cost function is non-convex. A sampling-based approach 
was proposed to capture the global minimum of the proposed cost function. Simulation results were also presented for both 
scenarios with the single-target and multiple targets occupying the same resolution cell. The obtained results confirmed the 
superiority of the optimal configuration compared to the common ULA structure in both single and multiple target scenarios. 
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Figure 13. The calculated cost for the two-target scenario. The cost was calculated for different values of A8 and scenarios with the optimal and ULA 
structure. 
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Figure 14. Location RMSE and CRLB for both optimal and ULA structures. The results are obtained for scenarios with different number of targets being 
located inside the same resolution cell. 



Appendix A 
Proof of Proposition^] 



Consider the n-th target and the m-th target that are located in the C!-th cell and the C2-th cell, respectively. First, assume 
that ci = C2 — 1. The unknown term Jg^-g^ is first calculated where the proof for other terms is similar. The entry Jg^-g^ 
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can be found using the following equality: 



( dp* \ ' y -i ( dp. 



T 



(A.68) 



Based on the definition of the covariance matrix in dl3t . it is evident that the second term on the right-hand side of the above 
equation is zero. Now, according to the definition of in d34l . the following equations can be derived for the partial derivative 
terms in ( IA.68I 1: 



9p* 
962 
dp* 
962 



962 062 892 362 

n d^(fj Cl ) d%(f) ci ) d^(fj C2 ) d$(fj C2 



962, d6m 992 992 
After some algebraic operations, the FIM in ( IA.68I ) can be written in the following form: 

'd^(f}ct)\ ( d$t(fj Cl )\' 



(A.69) 
(A.70) 
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( 8mc-i) \ f dU(fj C2 ) \' 
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99% 
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(9m*y\ ( dm^ v ( d$(f] ci ) \ ( dz(fj ca ) \' 



962 ) V 992 1 1 1 ! ' ' iA ' /l} 

Using ( [Tol l and ( fTTT > and the notations given by definitions [7] and [8] the following expressions can be derived for the derivatives 
in ( IA7TT I: 



V d62 ) \ 



961 



£(p?) / n(:.0n'(:,0PS (e «*«(!)) - C sin«(Z))) (C 2 cos(a#(Z)) - C™ 2 sin« 2 (0)) 



. 1=1 



962 ) \ 962 



$>*)'n(:,l)n'(:,0p£ sin«(0) +Q cos«(0)) sin(< 2 (0) + C 2 cos(w£(i))) 



(A.72) 



Therefore, the first term on the right-hand side of ( 1A.7U can be written as 



( 9m ci ) \ ( dK(fj Cl ) \' ( 9%(fj Cl ) \ ( 8Q(fj Cl ) \' 



\ 962 ) V 96 



962 J V 96 c n 



K 



P2 (1 - P2) 5>» )'«(:, On'O. OPS «" 2 cos« (/) - w « (0) + C' 2 sm« (0 - (0)) f (A-73) 



Other terms on the right-hand side of ( IA.7 11 > can be similarly found. The final form can be now written as follows: 

2 ( nm ~\ 



K 



^(p^)'n(:,0n'(:,0P^ cos «(0 -<(0) + t™ sin«(Q -<(0))| x 

" + Ml ~ /?,?)(! - P2) + Wffi + k 6 (32(l - 132)) (A.74) 



. z=i 
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where the term in the second line of the above equation can be written as follows: 



fai%(i-ffi)+h(i-ffi)(i-ffi)+hffiffi+w(i-i3%))= 08?)' o s- 1 omy 



(A.75) 



which is the coefficient Cg"i g c 2 . For the case with c\ = Ci, the same procedure can be followed and the expression in the 
proposition is similarly found. 



Appendix B 
Proof of Theorem[T] 

We begin with the optimization formulation given by d57| i. Define the new matrix T nm and the new variable t nm with 
{to =!,■••, M}, {n = 1, • • • , N}, and rewrite the optimization problem as follows: 



maX {ASn,--,AS„ m ,t,T*} 

S.T 



|As nm ||2 — d nrn 
|As nm ||2 ^ &nm 



(B.76) 



£™=i»*n + X£=i»r» = o 

(p ) ^nmP ^ ^nm 

(As nm )(As nm )' h T nm , V m = {!,■■■ , M}, n = {!,-■■ , N} 



where t = [tn ■■■ t nm ]', and T* = {Tn, ■ ■ ■ , T nm }. The second-norm terms in the constraints can be written into the 
following form: 
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(B.77) 



In addition, using the Schur-complement of a square matrix, the last constraint in (IB. 76b is written as 



i As: 



As T 
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(B.78) 
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Inserting the new forms provided by flB.77t and flB.78t in ( IB.76I ) and using the fact that (p c )'T„ m p c 
P = p c (p c )', the following new form is derived for the optimization problem: 

m=l 2^n=l lnm 



T{T nm P) with 



maX {T*,S*,t} 

S.T 



— 12x2 

As' 



•2 



do, 



^2x2 As 

As' d 2 

J nm nrn 



-< o 



(B.79) 



1 As 



/ 

urn 



^0,Vm = {!,•••, M},n = {!,•••, N} 



As T 

where 5* = {sti, • • • , s r Ar}. Now, the difference vector is written as 

As nm — s rn 

Replacing the above equation in ( IB.79K the form given by the theorem is obtained. 



(B.80) 



Appendix C 
Proof of Lemma[3] 

Consider the optimal structure found for the case with 0\. Assuming 62 as the new DOA the cost function in d57| i can be 
rewritten as follows: 

M N 

^^A S ; m P*As nm (C.81) 

rn — 1 n — 1 

with P* = (p*)'p* and p* = [cos(0i + A9) - sin(0i + AO)]'. The vector p* can be expanded as 
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(C.82) 



cos(A6») sin(A6>) 
- sin(A0) cos(A6>) 

Defining the first term on the right-hand side of the above equation as Gag, the cost function can be rewritten as follows: 

M N 

2~2 ^' nm G Ae PG' Ae As nm (C.83) 

m—l n—1 

We know that As°^ = (s°^ — s°J>) maximizes the cost function in (l57l i where P is the matrix corresponding the target with 
9\ as the DOA. Therefore, an optimal solution of the optimization problem with 62 as the DOA of the target can be obtained 

as 
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G'AO s t 
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S °l Vm = {l,---,M},n = {l,---,iV} 



(C.84) 
(C.85) 
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Consequently, the new optimal solution is written as follows: 

S°m = G A6S°m (C.86) 

s° 2 n = GAes°n, V m = {1, ■ • • , M}, n = {1, - ■ ■ , N} (C.87) 
Now, we have to check whether the new solution holds in the constraints. It can be shown that: 

||As£,||2 = ||AC|| a (C88) 

In addition, it is known that: 

M N / M N \ 

m—1 n—1 \m — 1 n— 1/ 

which implies that the new optimal solution also meets the constraints. 

Appendix D 
Proof of Proposition 0] 

Consider the optimization problem in dSTT i without the constraint on the mass center. In this case, the cost function is quadratic 
with respect to the unknown difference vectors. The unique optimal solution obtained by solving the resulting optimization 
problem can be written as {As° m } with m = {1, • • • , M} and n = {1, • • • , N}. It is evident that there are an infinite number 
of location solutions for which the above set of difference vectors are obtained. Let us define the i-th and the j'-fh sets as 
{S° l , S° 1 } and {S° J , Sy 7 '}, respectively. It is known from the geometry that: 

s# = G e s°J + bt, (D.90) 

where Gg denotes a rotation matrix with 6 as the angle of the rotation, and b refers to an arbitrary translation. Note that the 
above equation can be written for every other antenna in the array of receivers as well. Considering the mass center constraint 
given by ([3), we show that the translation should be zero in dP.90t . To show this, we first assume that there is a nonzero 
translation as bt;. Then, it is observed that such an assumption leads to the contradiction. From the assumption [3] it is known 
that the center of the mass of the array is located in the origin. Therefore, there should be another translation btv where 
btv = — ha- Under the new translations, the new difference vector is written as 

= ("3* -s° t i + 2b a ) (D.91) 

It is now evident that the new configuration gives a different set of difference vectors, which is a contradiction to our initial 
assumption (e.g., the same set of difference vectors). Therefore, the translation part in dP.901 > is zero. Now, consider the rotation 
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part in ( 1D.90I I. It is known that the rotation transform does not change the distance between each two points. Rewrite the cost 
function in ( |57T > into the following form: 



N N N N 



E E( s °-)' pCAs «™ = E E«™)' G « pCG » As «™ ( D - 92 > 



m=l n—1 m— 1 n— 1 

Using formal matrix operations, the new matrix U = G' g P c Ge can be written in the following form: 



U : 



(D.93) 



cos 2 (6» + 6» c ) - sin(0 + 6> c ) cos(6» + 6 C ) 

- sm{e + e c ) cos(e + e c ) sm 2 (e + e c ) 

The cost produced by each of two sets of optimal solutions is equal if the following condition is held: 

M N 

^^(AC)'MC-O (D.94) 

m— 1 n—1 

with V = U — P c . The equality in JD.94I ) is valid if either V — or V is neither positive nor negative semi-definite. First, 
assume V = 0. Based on the given form in flD.93t for the matrix U, it can be inferred that U — P c when 8 = nir. In other 
words, a rotation with mr as the angle of rotation provides the same cost function. Now, assume the other case where V ^ 0. 
It can be shown that matrix V has two eigenvalues {A, —A} where the value of A depends on the rotation angle and 9 C . 
Therefore, the zero inequality in flD.94t leads to a number of solutions for the difference vectors. The rotated configuration 
can be then another solution of the optimization problem if {S° l ,S° 1 } belongs to the set of solutions of JD.941 >. The above 
discussions state that the optimization problem provides at least two solutions for the optimum configuration of antennas. 
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